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I. INTRODUCTION 



The Standard Model (SM) of high energy physics suffers from the gauge hierarchy problem 
and the flavor problem. The first is the fine tuning required to maintain a low electroweak 
mass scale (Mew) i n t ne theory, in the presence of a high scale, the Planck Scale (Mpi). 
The second problem is a lack of explanation of the mass hierarchy and mixings of the quarks 
and leptons. 

Supersymmetry (SUSY) eliminates the gauge hierarchy problem by introducing for each 
SM particle, a new particle with the same mass but different spin. For example, for each SM 
quark/lepton a new scalar (squark/slepton), and for each SM gauge boson a new fermion 
(gaugino), is introduced. If SUSY is realized in nature, the fact that we do not see such 
new particles, we believe, could be because SUSY is spontaneously broken, making the 
superpartners heavier than the mass ranges probed by experiments. Owing to a lack of 
understanding of how exactly SUSY is broken, a phenomenologically general Lagrangian, 
for example, the Minimal Supersymmetric Standard Model (MSSM), is usually considered 
to compare with data. Various experimental searches have placed constraints on the masses 
and couplings in the MSSM. 

Attempts have been made to address the flavor problem by proposing various flavor 
symmetries. In a supersymmetric theory, a flavor symmetry in the quark sector might imply 
a certain structure in the scalar sector, leading to definite predictions for flavor changing 
neutral current (FCNC) processes on which experiments have placed severe constraints. 
In the literature, a lot of attention has been devoted toward analyzing the minimal flavor 
violation (MFV) scenario, in which the scalar flavor structure is aligned with the quark 
sector so that the two are simultaneously diagonalized. In MFV, the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix describes the flavor changing interactions in the supersymmetric 
sector as well, and the only CP violating phase is the one in the CKM matrix. In this work, 
we do not assume such an alignment, and we will consider non- minimal flavor violation 
(NMFV), which we treat as a perturbation over the MFV case. 

In this paper we wish to explore in what form a supersymmetric extension of the SM, 
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with a U(2) flavor symmetry, could influence K and B physics observables. We thus 
restrict ourselves to the quark and scalar-quark (squark) sectors. We consider an "effective 
super symmetry" []J framework, with heavy (TeV scale) first two generation squarks, in order 
to escape neutron electric dipole moment (EDM) constraints. This allows the possibility 
of having large CP violating phases in the squark sector. We consider a supersymmetric 
U(2) theory [2, E], impose recent K and B meson experimental constraints and obtain 
predictions for Bd — > X s j, Bd —>■ X s g, Bd — > X s £ + £~, Bd —>■ <f)K s , B S B S mixing and the 
dilepton asymmetry in B s . Though we consider a specific flavor symmetry, namely U(2), our 
conclusions would hold for any model with a sizable off-diagonal 32 element in the squark 
mass matrix. 

Some B physics consequences in a supersymmetric U(2) theory have been considered in 
Ref. P). Large tan/3 effects in B decays have been carefully analyzed in Ref. j^J, but for 
simplicity we will restrict ourselves to the case when tan (3 is not too large. Other work along 
similar lines, though in more general contexts, have been presented in Refs. jil If], B, Isj] . In 
this work we will include all dominant contributions to a particular observable in order to 
include interference effects between various diagrams. This has not always been done in the 
literature. We will then study the implications of recent data from the B-factories, including 
the b — > s penguin decay mode Bd — > <f)K s which shows a slight deviation from the SM 
prediction. 

The paper is organized as follows: In Section |H] we specify the supersymmetric U(2) 
theory we will work with, and the choices we make for the various SUSY and SUSY breaking 
parameters. In Sections IIHI and llVI we consider AS = 2 (Kaon mixing) and AB = 2 (BdBd 
and B S B S mixing) FCNC process, respectively. In Section|Swe will consider the implications 
of such a theory to AB = 1 FCNC processes, namely Bd — > Xsj, Bd — * X s g, Bd — > X s £ + £~ 
and Bd — > 4>K S . We conclude in Section IVII We give details of various squark mixings 
and their diagonalization in Appendix El and collect loop functions that we will need in 
Appendix El 
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II. SUPERSYMMETRIC U(2) 



The Model 



The supersymmetric model that we will discuss is as described in Ref. jjj, with the first 
and second generation superfields (ip a , a=l,2) transforming as a U(2) doublet while the third 
generation superfield (ip) is a singlet. The most general superpotential can be written as *: 



tab 



W = ipaxHip + —ipa 2 Hip a + —ip a a 3 Hilj b + 



S' 



ttb 



ip a a 4 Hip b + ^i! a a h Hi! b + /j,H u H d , (V 



where M is the cutoff scale below which such an effective description is valid, the «j are 
0(1) constants, and three new U(2) tensor fields are introduced: cj) a a U(2) doublet, <p ab a 
second rank antisymmetric U(2) tensor and S ab a second rank symmetric U(2_) tensor. The 
parameter [i could be complex and we allow for this possibility. Following Ref. 15 |, we assume 
that U(2) is broken spontaneously by the vacuum expectation value (VEV)t [3] 

(0 a > = f ° ) ; (0 a "> = ^ ab ; (S 11 ' 12 ' 21 ) = 0, (S 22 ) =V, (2) 

with V/M = e ~ 0.02 and v/M = e' ~ 0.004, in order to get the correct quark masses. 
These VEV's lead to the quark mass matrix given by (we show only the down quark mass 
matrix after the SU(2)l is broken by the usual Higgs mechanism) 

(d L \ 

C D -{d R s r R b R ) M d 



sl 



+ h.c, 



(3) 



In the superpotential each term encodes the "vertical" gauge symmetry, which, at the weak scale, is 
SU(3) x SU{2) x C/(l). Thus labels generations), 

ipiaHifjj = a u QiUjH u + a' a UiQjH u - adQiDjH d - a' d D?QjHd + (Lepton sector). 



t The dynamical means by which this VEV is generated is left unspecified. In general, (5 22 ) can be different 
from (4> a ), but for simplicity we will assume that they are the same. Also for simplicity, we take e, e' to 
be real. 
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where Vd = (hd) is the VEV of the Higgs field. In M. d) the Aj's are 0(1) (complex) coefficients, 
given in terms of the ati's. Ref. shows that such a pattern of the mass matrix explains 
the quark masses and CKM elements. 



If U(2) is still a good symmetry at the SUSY breaking scale, and broken (spontaneously) 
only below the SUSY breaking scale, the SUSY breaking terms would have a structure 
dictated by U(2). For our purposes it is sufficient to consider the down sector squark mass 
matrices, and they are given as 

(d L \ 



CD-(d* L s* L b\ 



. d*R S *R 



b R )M 2 RR 



( 



+ 



d* R ~s* R b* R )M RL 



fd L \ 



(d R \ 
\ 

h.c. 

/ 



(4) 



M\ L = M d M d + 



-ie'ml 



lemi 



V o 



M 2 RR 



M 2 RL 



MdMl 



7717 



• / 2 

-ie mi 



■ I 2 

ie 777,5 



V 



/i* tan [3M.d + Vd 



( O 
A ie > 
V o 



A 2 e 
A' 4 e 



\ 



2* 



em 



m 2 J 



+ D — term, 



LL 



\ 

+ e z m2 em 4 



+ D — term, 



ml ) 

\ 

A 4 e 

A3 J 



RR 



(5) 



where m\ and A4 are determined by the SUSY breaking mechanism. Here m\ , m|, m\ and 7775 
are real, while 7774 and Ai could be complex. We will assume that the Ai are of order A, a 
common mass scale. The D terms are flavor diagonal, and since we are interested in FCNC 
processes, we will not write them in detail, but will think of them as included in m\ and 7773. 
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Thus far we have presented the mass matrices in the gauge basis. In the following sections, 
we will work in the superKM basis in which the quark mass matrix is diagonal, and the quark 
field rotations that diagonalize the quark mass matrix are applied to the squarks, whose mass 
matrix would also have been diagonalized in the MFV scheme. Since we will not assume an 
MFV structure, in the superKM basis, there would be small off-diagonal terms in the squark 
mass matrix, which we treat as perturbations. The structure of the squark mass matrix in 
the superKM basis is similar to that in Eq. owing to the smallness of the mixing angles 
that diagonalize the quark mass matrix. 

B. SUSY parameters 

Lacking specific knowledge about the SUSY breaking mechanism realized in nature, 
we make some assumptions on the SUSY mass spectrum. Neutron EDM places strong 
constraints on the CP violating phases and the masses of the first two generations of scalars. 
To satisfy this and other collider constraints, we consider an "effective SUSY" framework in 
which the scalars of the first two generations are heavy, suppressing EDM, and allowing for 
larger CP violating phases. Defining the scalar mass scale, mo ~ 1 TeV, we take all rrii ~ uiq 
except for i = m 3RR ~ 100 GeV. We take A ~ m , the gaugino mass parameter M 2 
and charged-Higgs masses to be 250 GeV and the gluino mass to be 300 GeV*. We assume 
such a spectrum just above the weak scale without specifying what mechanism of SUSY 
breaking and mediation might actually give rise to it. As we will show later, if realized in 
nature such a spectrum would lead to enhancements in the processes we are considering here. 

The rates of various FCNC processes follow from the mass matrix that we have specified 
in Eq. (jSJ). We will work in the superKM basis. The interaction vertices in the mass basis are 
obtained by diagonalizing the mass matrices in Eq. (j5J), and the perturbative diagonalization 
to leading order is shown in Appendix. [X] 

* The Tevatron bounds on the stop, sbottom and gluino masses are discussed in Ref. We note here that 
the bounds in general get less stringent as the neutralino mass increases. 
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The dominant NMFV SUSY contributions to FCNC processes would be due to the 32 
and 23 entries in Eq. since they are the biggest off-diagonal terms. For convenience we 
define 

cRL,RR,LL (•^RL.RR.Ll) '32,23 ,. 

0? 9 ^ = ! — k • (o) 



'32,23 



Since we have written down an effective theory and not specified the dynamics of U(2) and 
SUSY breaking, we can only specify the order of magnitude of 532,23- To parametrize this 
uncertainty we write, 



s:RL _ v dAe RL 
°32,23 ~~ - 1 u 



m, 



2 "32,23' 




cLL,RR _ em 4 jLL,RR 
"32,23 ~~ HI fl 32,23 > 



mo 



(7) 



where we have denoted the unknown 0(1) coefficients by d 32 



LL,RR,RL 



23 



m 


1000 


tan f3 


5 


OR,tR 


100 




200 e i2 - 2 


d.R,SR 


1000 


M 2 


250 


mq L 


1000 


M~ g 


300 


A 


1000 


m H ± 


250 




2e 43 - 2 


dW 


1.75 e iL6 



TABLE I: Default SUSY parameters for this work that satisfy all experimental constraints discussed 
in this paper. All masses are in GeV. 

We summarize our choice of the parameters in Table |1J For these values, from Eq. (|7J), the 
natural sizes of 5 32 23 are given by 



<e 23 = 6.82 xl0- 4 dg 23 , 



^23^ = 0.02 4^ 



We will find in the rest of this paper that 5 RL induces NMFV AB = 1 FCNC processes 
dominantly, while S RR,LL induces AS* = 2 and AB = 2 FCNC processes. Though the 5 R 2 L 
and 823 elements have similar magnitudes, the gluino NMFV contribution to AB = 1 
FCNC processes is larger, since we take &r to be much lighter than the other scalars, and 
the gluino diagrams are relatively suppressed by the heavier hi mass. Therefore, in 
this work we will include only the dominant S R 2 L contribution. We illustrate this in Fig. ^ 




FIG. 1: Gluino contribution to — > X s ^. The diagram on the left, proportional to has the 
lighter scalar bji, while the one on the right, proportional to 8^, only has heavier scalars and is 
therefore relatively suppressed. 

where we show the gluino contribution to Bd — > Xgj as an example. Similarly, owing to the 
smaller &r mass, the S§ 2 ^3 NMFV contribution to AS = 2 and AB = 2 FCNC processes is 
relatively larger compared to the S^ 2 3 contribution. We note here that, from Eq. (jAll|) in 
Appendix El the sbottom mixing angle is negligibly small, and therefore, we ignore sbottom 
mixing effects; stop mixing is not as small and we include its effects. 

In the next three sections we will discuss the implication of the U(2) model to AS* = 2, 
AB = 2 and AB = 1 FCNC processes. From this we will see that present experimental 
data are compatible with the values shown in Table HI and we will obtain expectations for 
some measurements that are forthcoming. We will present plots of different FCNC effects by 
varying a couple of parameters at a time, while keeping all others fixed at the values shown 
in Table 



III. AS = 2 FCNC PROCESS 

The CP violation parameter due to mixing in the Kaon sector has been measured to 
be Q 

\e K \ = (2.284 ± 0.014) x 10~ 3 . (9) 

We wish to estimate the new physics contributions to ex in the scenario that we are 
considering. Here we note that even though the direct CP violation parameter e' K /eK has 
also been measured, large hadronic uncertainties do not permit us to constrain new physics 
models through this observable. 
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Kaon mixing is governed by the AS = 2 effective Hamiltonian 

5 3 



^ = 1^ + 1^. ( 10 ) 

i=i i=i 

where, 

Q3 = ^H S L^R S L 5 

Q4 = d^sf^d^s^ , 

Q 5 = dS^aS • (11) 

The operators Qi (i=l,2,3) are obtained by exchanging L <-> i?. In the SM and the new 
physics model we are considering, the dominant contributions are to Qi, as we explain later 



in this section. The CP violation parameter ex is then given by (see for example Ref. 121 ]) 



e K = e / —-^—^[C^)) , (12) 

where Bk is the Bag parameter and fx is the Kaon decay constant. 

In addition to the SM W box diagram contribution to C±, in the supersymmetric U(2) 
theory we are considering, the charged-Higgs and chargino MFV contributions could be 
sizable. The dominant MFV contributions to C\ can be written as 

Cf FV = Cf + Cf + Cf , (13) 

which is the sum of the SM W, the charged-Higgs, and the chargino contributions, 
respectively. 



SM contribution: The SM W contribution is 11] 



C! M {m t ) = CY{m t ) = ^L |(^ s ) 2 S (x t ) 

+ (v; d v cs ) 2 s (x c ) + 2 (v td v ts v; d v cs ) s (x u Xe ) ) , (14) 
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where the function Sq is given in Appendix iBl Eq. (|B3|) . and x t = mf/m^y, x c = m^jm^. 
The QCD correction due to renormalization group running from m t to mj, gives 

Cf M {m K ) = CY(m K ) = {(Oy Vsa S (x t ) 

+ (V: d V cs ) 2 VK22 S (x c ) + 2(V t * d V ts V; d V cs )r ]K32 S (x t ,x c )\ , (15) 

where the i]k are QCD correction factors given in Eq. (J2*U|) below, and are the CKM 
matrix elements. 

Charged-Higgs contribution: Supersymmetric theories require two Higgs doublets to give 
masses to the up and down type fermions. The Higgs doublets contain the charged-Higgs 
H ± , and the dominant charged-Higgs-top contribution is jl^l p 

C*{m) = ^ (KM 2 [-F*] , (16) 



1 1 2 

F v = 77 — TB x t Y ii r H, r H , x t , x t ) + — — y-zx\ Y"i(l, r H , x t , x t ) - - — ^xt Y 2 (l, r H , x t , x t ) 
4 tan 4 p 2 tarn p tan z p 



where r# = rn^/m^r, and the functions Y\ and Y 2 are given in Appendix [TJ Eq. ijEl) . 
Chargino contribution: The dominant chargino-right-handed-stop contribution is \l 



pX _ llpW |2lpU; |2 y,/ \ 



where rj 2 = m| /mfp, si^ = 2 /myy, and the coupling is given by 



r(0 -./o(r*\ <r*\ !^kEk m ' ris^i 

r xB - V2(^) u (C f ) 12 — ^— — . ( 18 ) 



with the chargino and stop diagonalization matrices (Cr) and (Cf) given in Appendix \K\ 
Eqs. QA6[) and (|A10|) . respectively . Taking into account renormalization group running, we 
have 

C^(m K ) « to3 < » X K) . (19) 

Gluino contribution: In general, the NMFV gluino contributions induce many operators 
shown in Eq. (JUJ), but in the model we are considering, these are not significant due to a 



The charged-Higgs also contributes to the operator Q2, which becomes important only at large tan/?. 
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suppression from the heavy d and s masses, Glashow-Iliopoulos-Maiani (GIM) suppression 
owing to their approximate degeneracy (split only by C(e 2 ), cf. Eq. ©), and the contribution 
from the relatively light right-handed sbottom being suppressed by its small mixing to the 
first two generations. Moreover, owing to the structure of the mass matrix, Eq. (J3J, the 
gluino contribution is real, and hence does not contribute to ex- 

In our numerical analysis, we take the following values for the various parameters jlfl. 11]: 



VK33 = 0.57 , T] K2 2 = 1.38 , T] K32 = 0.47 , 

f K = 0.160 GeV , 0.6 < B K < 0.9 , (20) 
m K = 0.497 GeV , Am K = (3.48 ± 0.01) x 10~ 15 GeV , m c = (1.2 ± 0.2) GeV . 

The SM prediction for ex is in agreement with the experimental data, but it should be 
noted that there is considerable uncertainty in the lattice computation of the Bag parameter 
B K (see Eq. (|20j)). The chargino and charged-Higgs contributions to C\ add constructively 
with the SM contribution. Therefore, if the true value of Bk is taken to be closer to the 
lower limit, we can allow MFV contributions to be up by a factor of 1.2 compared to the 
SM value; i.e., lm(Cf IFV ) /lm(Cf M ) < 1.2. Fig.Elshows the region of MFV parameter space 
where this is satisfied. This justifies some of the choices we make in the list shown in Table |U 



IV. AB = 2 FCNC PROCESSES 



A. General formalism 



We start by d ls cu Ssillg in genera! B,B, nnxing and later s pecM 1Z e in succession to B d B„ 
(q=d) and to B S B S (q=s). The AB = 2 effective Hamiltonian is given by m- 

n e /L2 = Y, c iQi+ y L <5 iQi > ( 21 ) 

i=l i=l 

where, for B q , 

11 




FIG. 2: The (dash-dot, dash, solid) curves are (1.05, 1.1, 1.2) contours of Im(C 1 MFl/ )/Im(Cf M ), 
showing the MFV contributions to Kaon mixing relative to the SM. Parameters not shown on a 
plot's axes are fixed as shown in Table 



Q 2 


= m&L , 
















Q 5 




(22) 



The operators Qi (i=l,2,3) are obtained by exchanging L ^ R. The Wilson coefficients Cj 
are run down from the SUSY scale, M s , using [3] 

C r (m b ) = E 0^ + V^) ifC.(Ms) (23) 

i s 

where r\ = a s (Ms)/a s (m t ) and the cij, fej and c, are constants given in Ref. 



The matrix elements of the Qi in the vacuum insertion approximation are given by |l3l . Il5 ] . 

2 



2 



fl.IO.WIB.) = 4 P^-] 



"9/ 
2 



1 / m B„ \ 2 ,2 



12 



B q \QM\B q 
B q \QM\B q 



2 \m b + m q 



m B fB a B ^) 



1 / rn Ba 



™B q fl q B M 



(24) 



6 \mb + vn q t 

where we take for the decay constants f Bq = 0.2 ± 0.03 GeV and the Bag parameters (at 
scale m b ) B x = 0.87, B 2 = 0.82, B 3 = 1.02, B 4 = 1.16 and B 5 = 1.91 HQ- 

The B q mass difference is given by 

Am Bq = 2\M 12 (B q )\ , (25) 
where Mi 2 (B q ) is the off-diagonal Hamiltonian element for the B q B q system, and is given by 

Mia = Mf 2 M + M* 2 USY , 



r 



12 ~ r 12 



SM 



Sto an excellent approximation is dominated by the SM tree decay modes. From Refs. 
we have, 

1 



iSM 
12 



\M 12 (B q ) 
G 2 F m\m Bq B B Jl q 



2m Ba 



BJTt 



8tt 



AB=2 




1 2 

4 


1 


4 3 

-zi - 

3 c 





(26) 



where v x = V xb V* z c = m c /mi and we take £?b 9 « 1.37. 
The dilepton asymmetry in B q is given by Q| 

. o _ iV(B,B,) - N(B q B q ) 



A q = 



N(B q B q ) + N(B q B q ) 



Im 



/r 



12 



VM 



12 



(27) 



We discuss next the SM and new physics contributions to the coefficients Cj and Cj. 
MFV contribution: The SM contribution is almost identical to that shown in Eq. ((Hj) 
but for the fact that it is sufficient to keep only the top contribution (the So(%t) term) and 
changing the CKM factor to (v^V^) . The new physics MFV charged-Higgs and chargino 
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contributions are again identical to Eqs. (fTBj) and (fTTj) . respectively, with the same change 
for the CKM factors. Ci(m t ) is evolved down to m& using Eq. 



Gluino contribution: We only include the dominant gluino-right-handed-sbottom box 
diagrams with 8§f and 8§P mass insertions, since bn is the only relatively light down type 
squark in our scenario. These contributions are given by 



Cf(M- g ) 
Cl(M s ) 

c!(m s ) 



Ws 



Ws 



r RR f r RR\ * 
1 63 \ L q3 ) 

r RL f r RL\*' 
1 63 \ L q3 ) 

r RL fr^ L \ 
1 63 \ L q3 ) 



— I A + -M?I 4 

.36 9 9 

2 9 



1 

-M n 
2 9 



M?I 4 



with the box integrals I 4 and I4 given in Appendix [5] The couplings are given by 



r RR 
1 63 



cos#f 2 R 



T* R = sin 0* R sin e^e-*^ - ' - 



= sin^ 9 L sin#£V 



'12 



7 32 



rf 3 i = cos^t 

1 s3 — 



cos^sin^e"^ 
rf 3 i = -cosCsin^ i e-^. 



(28) 



(29) 



obtained from the 3x3 mixing matrix that is the product of C^ r - r and Ci rSr , with the 
mixing angles 9 and phases 7 given in Appendix El In our U(2) model, if m 4 is of the same 
order as A, based on the estimate in Eq. (|SJ), we expect Cf to receive the dominant gluino 
contribution from 8§f'. We will focus on this contribution in the following. 



We point out in Appendix |XJ Eq. (|A18|) . that (IrSr mixing can be generically large (near 
maximal), in which case the gluino contributions to both B^Bd and B S B S mixing can be 
sizable. However, if e'm\ -C £ 2 m\, this mixing can be small and the gluino contribution to 
BdBd mixing is negligible since it is proportional to sin^f^, cf. Eqs. (I28j) and ()29)) . The 
gluino contribution to B S B S , however, can still be sizable in either case since it is proportional 
to cos 6? 2 R . 
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FIG. 3: The (dash-dot, dash, solid) curves are (1.1, 1.2, 1.3) contours of \C^ IFV /Cf M \, showing 
the MFV contributions to B^Bd mixing relative to the SM. Parameters not shown on a plot's axes 
are fixed as shown in Table [fl 

B. B d B d mixing 



The BdBd mass difference (Amj, and CP violation in Bd — > ipK s (a^J have been 



measured to be 
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Am d = 0.502 ± 0.007 ps" 1 , 

a^ Ka = 0.725 ± 0.037 . (30) 
In the SM, the usual notation is, = sin 2/3. 

As we have already pointed out in Section IIII[ the charged-Higgs and chargino MFV 
contributions add constructively with the SM contribution. The SM prediction agrees quite 
well with the data, but given the uncertainty in fg d , cf. below Eq. (124)) . it might be possible 
to accommodate an MFV contribution up to a factor of about 1.3 bigger than the SM 
contribution. We show in Fig. |3] the region in MFV parameter space that satisfies this 
constraint, ignoring the gluino contribution. 

As pointed out in the previous subsection, in general we expect in the U(2) model, gJrSr 
mixing to be near maximal, in which case the gluino contribution to BdBd can be sizable. 

15 



The gluino contribution can then be important to both Arrid and a^K ■ Taking this into 

I I 

account, we can write cl^k, = sin {2f3 + 26 d ), where 6 d is the new phase in Mi 2 (-E>d) 18], and 
we have 1G, 1911 



a>ipK s — Im(A^) , 




q A(B d ->iPK s ) 
P A(B d -^^K s ) 



v \M W - f r 12 ' 



(31) 



with M 12 and T 12 given in Eq. (|26p. (The "— " sign in X^k is because the final state is CP 
odd.) In our case, T 12 <C M 12 , so that 



where "arg" denotes the argument of the complex quantity. 

For the case when c^rSr mixing is large, we show the gluino contribution to B d B d in Fig0] 
The plot on the left also shows the constraint from a^K s , which is not shown in the plot on 
the right since almost the whole region shown is allowed. The region (rr < arg (5§f") < 2n) 
is not shown since it is identical to the region (0, n). From the figure, we see that in the 
large mixing case, the constraint on 5^ is quite strong. However, if (IrSr mixing is small, 
the constraint on 5§P from B d B d mixing is weak. 

C. B S B S mixing 

B S B S mixing has not yet been observed and the current experimental limit is Am^ s > 



14.4 ps~ l @ 95% C.L. [10j. The SM prediction is: 14 ps" 1 < Am Bs < 20 ps" 1 |2fl|. The SM 



prediction for the dilepton asymmetry Af t s is small, around 10 4 , cf. references in Ref. [16]. 

B S B S mixing depends quite sensitively on 6§p, and for the region in Fig. 0] allowed by 
B d B d mixing, we find Ame s 22 ps _1 and Af^ « 5 x 10~ 4 . This Am^ is a little higher 
than the SM prediction, although may be within the SM allowed range, given uncertainties. 

As we pointed out in the previous subsection, if (IrSr mixing is small, then the B d B d 
mixing constraints on becomes weak. If such is the case, there are essentially no 



a^K s ~ sin(arg(Mi 2 )) , 



(32) 
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FIG. 4: For large (Irsr mixing, the (dash-dot, dash, solid) curves are (0.9, 1.0, 1.25) contours of 
\{Ci + Ci)/Cf M \, showing the NMFV contributions to B d B d mixing relative to the SM. d§^ is 
defined in Eq. ((7J). The hatched region is excluded by cl^k s - Parameters not shown on a plot's axes 
are fixed as shown in Table [IJ 

constraints from B d B d mixing, and we show contours of Am# s and Af t a in Fig. |S| We 
show only the range (0 < arg (S§^) < it), since the (n, 2ir) range is identical to this. It can 
be seen that Am# s can increase significantly above the SM prediction. The projected Run 
II sensitivity for Am^ at the Tevatron with 2 fb _1 is around 40 ps _1 j^j], and can probe 
a significant region of U(2) parameter space. If a higher value of Am^ is measured than 
what the SM predicts, it would indicate the presence of new physics. Measuring A^ 3 can 
also significantly constrain 5§f~ as can be seen from Fig. El (right). 
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FIG. 5: For small (IrSr mixing, the (dash-dot, dash, solid) curves are (15, 25, 40 ps 1 ) contours 
of Am Bs (left), and (1(T 4 , 1CT 3 and 1CT 2 ) contours of |A^ S | (right). d§^ is defined in Eq. Q). 
Parameters not shown on a plot's axes are fixed as shown in Table [I] 

V. AB = 1 FCNC PROCESSES 
A. Effective Hamiltonian 



The AB = 1 effective Hamiltonian at a scale /i in the operator produce expansion (OPE) 

is 



G 



niLx = —kvtsv; b £ c^o^) + cv 7 (/i)o 7 » + c 8g ^)o 8g ^) 



with 



y/2 



(33) 



K i=l... 6,9,10 

01 = (s a C /3 )y_ J 4(c /3 6 a )y_A, 

2 = (sc)v-.a(c&)v-a, 

3 = (sb) V - A ^2(qq)v-A, 

1 

04 = {s a bp) V _ A J2(qp<la)v-A, 

1 
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5 = (sb) v ^ A J2( ( l ( l)v+^ 

g 

6 = {s a b(3)v-A^2{qpq a )v+A, 



q 

g 

°7 7 = ^^m h s a a^ u {l + 75)60-^/11/, 
9 = (s6) y _ A (ee)y, 

Oio = (sb) v _ A (ee) A , (34) 

where, the subscript (V ± A) means 7^(1 ± 75), and F^ u , are the electromagnetic and 
color field strengths, respectively. 

The Wilson coefficients can be computed at the scale M w (the W boson mass), and then 
run down to the scale m& (the b quark mass). Below, when no scale is specified for the 
coefficients, it is understood to be at m&, i.e., Ci = CArrih). The coefficients when run down 
from M\y to mix under renormalization, so that I23J ^ 

Cj = EM°'- (j = 1,...,6), 
i=i 

C 7l = ^C 7 ^M W ) + I (77M _ 77!) C 8g (M w ) + X; h tV a *C 2 (M w ), 



3 

C 89 = V ^C 8g (M w ) + £ h tV a >C 2 (M w ), (35) 



s 

14 



where 77 = w 0.56 and /ij, /ij, a, and kji are given in Ref. In addition, the 

evolution equation for Cg is given in Ref. j^j], and Ciq is not renormalized. 

Separating out the new physics contribution to the renormalization group evolution, i.e., 
Eq. ((US)), we get 

ri r*SM 

C 7l = + 0.67C^ W (M W ) + 0.09C^ W (M W ), 

Cs 9 = Cg g M + 0.70C™(M w ), (36) 



^ Here, as a first step, we use the leading order result. The next to leading order result can be found in 
Ref. Q. 
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in which the superscript "SM" indicates the contribution from the SM, and "new" from new 
physics. 



SM contribution: The SM contribution to C 7l (My/) and Cg g (Mw) are given by 
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26] 



qSM 



(M, 



w, 



w , 



-F 



LL 
7,8 



m t 



(37) 



(3? 



where F^(x) are given in Appendix |H1 Using Eq. (J35|) we can compute Qf M , C^ 1 and 
r<SM 

In the following, we will discuss, in order, the new physics contribution arising from the 
charged-Higgs boson (H^), charginos (x^) an d gluinos (g). 

Charged Higgs (H ) contribution: The charged-Higgs contribution to Bd — > X s 7 is 

given by 0, T 



26, 2, 



(M, 



1 



cot 2 (3F 7 L j 



LL 



(39) 



where F^(x) and Ff%(x) are given in Appendix IB1 

Chargino (x^) contribution: The chargino-stop contribution can be comparable to the 
SM contribution for a light stop and chargino. In the scenario that we are considering, the 
stop mixing angle is negligibly small and m; ~ rrif ~ mo and mr ~ m- t ~ My/- We 



therefore run the t\ contribution from % down to My/ and evaluate the t?. contribution at 



M w . The chargino-stop contribution is 5j, 
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3=1 
3=1 



I ^7,8 I 
A/f 2 

l F L I ^7,8 



ij M w pRL 

r 7,8 



M 



Xj 



m 



2 • M w M 



m 



fa 



(40) 



where the loop functions F^ are given in Appendix El and T l [ and Yrl contain the stop and 
chargino mixing matrices. Explicit expressions for T^, Yrl an d the renormalization group 
equations to evolve ^^(mo) down to M w are given in Ref. 
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Gluino (g) contribution: In our NMFV scenario, the gluino contributions can be sizable 
since they couple with strong interaction strength. Furthermore, because the sbottom mixing 
angle is negligibly small, m~ bL pa m~ bi ~ m n and mj. ~ rn 
enhanced piece, the gluino contribution is 



'62 



M w . Keeping only the — 



CL(M W ) 



-— — — 9 - cos 6>" L sin 6^e n ™ - 

G F V t * s V tb m b 32 32 9 

47TasV * **1 cos e« L sin 0§?e-+& I 
G F V t * s V tb m b 32 32 8 




(41) 



where the mixing angle 9§^ and phase 7^ L are defined in Appendix El and _F 4 and F g are 
defined in Appendix El In the above equation, we have neglected the effect of running the 
b\ contribution from rfiQ to Mw as the 62 contribution is dominant. 

The dominant new physics contribution is given by adding Eqs. (J39|) . (J4Uj) and ([41)1 . which 
yields 



(M, 



W , 



77,89 



(M, 



77,89 



(M, 



w ) 



(42) 



In what follows we will discuss in detail the new physics contribution predicted by the U(2) 
model to the rare decay processes B d — > X s 7, B d — > X s g, B d — > X s £ + £~ and B d — > 4>K S . 



B. B rt 



Xsl, B d ^ X s g 



The dominant operators contributing to B d —>■ Xgj and B d —>■ X s g are O2, and 
0$ g . The decay branching ratio B.K.(B d — > X.^), at leading order, normalized to the semi- 
leptonic B.R.(B d -> X c ez>) pa 10.5%, is given by ^bab^ 



T(B d -> X s7 ) 



r(5 d 



X c ez/) 



6a 



TT ^ 

3 v m;, ' 



^t6 



cb 



\r 1 2 

IW7I ) 



(43) 



(e 1 >(i-&)e™°-*) 

where g(z) = 1 — 8z 2 + 8z 6 — z s — 24z 4 ln(z) is a phase space function, and S is the fractional 
energy cut, i.e., only photon energy E 1 > (1 — 5)E™ ax is accepted. 



21 



The CP asymmetry in B d — > X s ^ is given by 

T(B d -> X sl ) - T(B d - X s - 7 ) 



A B d ^X a ~/ 



T(B d 
1 

|r"< 12 
I W7I 



x s7 ) + r(5, -> xrr) 



E y >(l-6)E» 



{a 27 (5)Im[C 2 C*J + a 87 (5)Im[C 8g C* 7 J + a 28 (5)lm[C 2 C* 8g ]} . (44) 



For 5 = 0.15, which is a typical experimental cut, we use 027 = 0.0124, 037 
a 2S = 0.0004 



-0.0952 and 



+0.3 \ 

-0.28 ; 



The experimentally measured |17J branching ratio is B.R.(B d — ► X s 7) = (3.52 
10" 4 . In the SM, we have C| M « 1.11, Cf 7 M « -0.31 and Cf/ 7 « -0.15. The SM 
prediction for B.R.(B d — > X/y), which depends on |CV 7 |, cf. Eq. fjlHj) . is largely consistent 
with experiment, and new physics contributions to |CV 7 | is constrained by this branching 
ratio. In the context of SUSY this has been analyzed, for example, in Refs. j^, 13, 0] • 

The SM CP asymmetry in B d — > X f 7 is of the order of 1%, so that a larger CP asymmetr 



measured would imply new physics 
-0.07 < Ac d j^ Xs ~< < 0.07. 



29j. The present limit at 95 % C.L. is [17, 
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The B.K.(B d — > is obtained simply from Eq. 



T(B d ^X s g) 
T(B d -> X c ei/) 



C(i2) 



7T 



?7lb • 



V, 



cb 



ic« 



(45) 



where the SU(3) quadratic Casimir C(R) = 4/3. The B.K.(B d — > has large 

experimental and theoretical uncertainties and Ref. suggests that the data might prefer 
a B.R. value of around 10%. 

Figs. EH and [3 show the interplay between the W^, H ± , an d g contributions to 
B d — > Xgj, where the sum of these contributions to the magnitude of C 7l is constrained 
by B.K.(B d —>■ X s j). The experimental data on B.K.(B d — > X^) allows (at 2a) the region 
bounded by the contours shown in the figures. In the plots, the parameters are as given in 
Table |1J and some relevant ones are varied as shown in the figures. 



To illustrate the dependence on the MFV parameters, we consider for example, in Fig. El 

,, — ► Xg'-f) as a function of tuh and M 2 (left) and 



the dependence of A^fp* Xa ~ i and B.R. {B, 



as a function of tan/3 and arg (/i) (right), for the choice of parameters shown in Table H] 
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FIG. 6: The boundaries between the shaded regions show (2.5,3,3.5)% (darkest to lightest) contours 



l CP 
l CP 

contours of B.R.(I?d — ► X s j). Parameters not shown on a plot's axes are fixed as shown in Tabled 



of A J ^ d p* Aal as a function of tuh and M2 (left), and, (-3.5,0,3.5)% (darkest to lightest) contours of 
Aqp^^^ as a function of tan j3 and arg (ii) (right). Superimposed is the experimental 2 a allowed 



The experimental 2 a allowed contours of B.R.(i?^ — > X^) are also shown. In Fig. (left), 
the shaded regions show A^ d p* Xsl as a function of the magnitude and argument of d§[, 
the dimensionless 0(1) coefficient defined in Eq (JSJ). In Fig. [7| (right), we show contours of 
B.R.(5^ — > X s g), and a B.R. of up to about 15% can be accommodated in this model. 

C. B d - X s ^+^- 



The dominant operators contributing to B d — > X s £ + £ (£ = e, /i) are O^, 9 and Oi . It 
is usual to define 



and 



cm 

Cio 
s = 



Z7T 



— ( 

2tt 



— wo j 



mi 



(46) 
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FIG. 7: The boundaries between the shaded regions show (-7,-3,3,7)% (darkest to lightest) contours 
of A^p* Xsn< (left) with experimentally allowed 2a contours of B.R.(l?d — > X s ^) superimposed, and, 
1 %, 7.5% and 15% contours of predicted B.R.(Z?d — > X s g) (right). Parameters not shown 
plot's axes are fixed as shown in Table [IJ 



on a 



The (differential 1 ) partial width -^T(B d 
by Q: 



X s £ + £ ), normalized to T(B d — > X c eu), is given 



r(5 d -> X c e^) 



a 

4^ 



cb 



(1-i) 



f( mc )rJ mc ' 



l + 2S)(|C 9 e// | 2 + |C 10 | 2 ) + 



4 ( 1 + - ) |C 77 | 2 + 12 Re(C 77 C' 9 e// ) 



.(47) 

where / and k are phase space functions and Cg is the QCD corrected Cg, given in terms 
of Cg and Cj (i=l...G) 22]. Integrating this we get the prediction for the decay branching 
ratios and we show this in Table |H] for the SM along with the experimental result H. We 
choose the lower limit on the integration to correspond to a typical experimental choice, 
(j>e+ + Pi-) 2 > (0-2 GeV) 2 . Since the rate of B d — > XJ + Pr is down by the square of 
the electromagnetic coupling constant compared to B d — > -^7, the experimental errors are 
comparatively larger. 
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Experiment [17] 


SM prediction 


B.R.(B d X a i+£~) 


4.46ig;|| x 1(T 6 


5.3 x 10" 6 



TABLE II: The current data for B d -> X s £ + £~. 



arg fd LR 32 ) 




FIG. 8: The (dash-dot, dash, solid) curves are (5.25, 6.25, 7.25) x 10~ 6 contours of B d -» X s £ + £-. 
Parameters not shown on the plot's axes are fixed as shown in Table [I] 

Fig. |H] shows the contours of B.K.(B d — > X s £ + £~) as a function of r. Compared to 
B.R,.(B d — > X s 7), cf. Fig. [7| (left), the B d — > X s £ + £ _ constraint is not very stringent 
right now, and improved statistics at the B-factories could place tighter constraints on the 
parameter space. 

D. B d -» cf>K s 

The decay B d — ► (6 — ► sss at the quark level) can be a sensitive probe of new 
physics since the leading order SM contribution is one-loop suppressed, and loop processes 
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involving heavy SUSY particles can contribute significantly. However, the computation of 
B.R,.(B d — > (j)K s ) suffers from significant theoretical uncertainties in calculating the hadronic 
matrix elements. We follow the factorization approach, details of which are presented in 
Ref. 35]. The theoretical uncertainties largely cancel in the CP asymmetry, and is therefore 
a good probe of new physics. 



The CP asymmetry in B d — > <pK s is defined by 

r (B d (t) - <pk s ) - r (B d (t) - <px s 



A 



CP 



r \B d {t) -> 4>k s ) + r (B d (t) -> <j>K s ) 

-C^ K cos (Am Bd t) + S^k sin (Am Bd t) , 



(48) 
(49) 



where 



S<f>K = 



A 



4>K\ 



l + |A x| 2 ' 
2lm(X 4>K ) 

1 + \^k\ 2 ' 

r -2i((3+e d ) A ( B d 



<pK s 



A(B d -> <f>K.) 

where B d (t) represents the state that is a B d at time t — 0, ArriB d is the B d B d mass 
difference, /? is the usual angle in the SM CKM unitarity triangle fits to the CP asymmetry 
in B d ^ J/ip K s , and 9 d is any new physics contributions to B d B d mixing (9 d is discussed 
in Section TlV B|) . The SM predicts that the CP asymmetry in B d — > <pK s and B d — > J/ip K s 
should be the same, i.e., Sqk — sin 2/3. 

The B d — > 4>K S amplitude and partial decay width are given by jsl l^: 



A{B d -> <PK S 



T{B d - 0AT s 



£ A 

p=u,c 

(■if:,'" : 



(a 3 + + as) — -(07 + a 9 + a 



32tt 



(F" ) \A{B d — > 0K S 



ml m 



K - 



2 • 



(50) 
(51) 



where the phase space function** X(x, y,z) = 
constant /<* = 237 MeV, the form factor Ff^ K 



x 2 + y 2 + z 2 — 2xy — 2yz — 2zx, the decay 



0.38 and A r 



VpbV* 8 . The SM Oi's, 



We thank Liantao Wang for clarifying the expression for X(x, y,z). 
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in terms of the Cj's, are given in Ref. 



. 3 



to which we add the new physics contribution 
given in Eq. (|42|) . We do not include the power-suppressed weak annihilation operators and 
we refer the reader to Refs. j3 an d 36 1 for a more complete discussion. As explained in 
Section. Ill ±31 we are only including the SUSY contribution, as this is the dominant one. 
The amplitude for the CP conjugate process B d — > (f)K s is obtained by taking \ p — > X* 



The current B d — > <pK s experimental average [17] is summarized in Table. IIHI The SM 
requires Sqk = Sj/^k = sin 2/5, but the experimental data has about a 2cr discrepancy 
between S^k an d Sj/^k-^ Though not convincing yet, this could be an indication of new 
physics and we ask if this can be naturally explained in the theory we are considering. 





Experiment [17] 


SM prediction 


B.R.(B d -> <f>K a ) 

S<f>K 
C(j)K 


8.3t\i x 10~ 6 

0.34 ±0.2 
-0.04 ±0.17 


~ 5 x 10~ 6 
0.725 ± 0.037 




TABLE III: The current data for B d -> <j>K a 



We showed in Section II V Bl that if (IrSr mixing is small, there is no significant new phase 
in Mi 2 (B d ) (i.e., 9 d ~ 0). For this case, we scan the parameter space arg(/i), \5§^\, arg(<5|| i ), 
and in Fig. El show a scatter-plot of the points that satisfy all experimental constraints 
including B.K.(B d — > X s j) and B.K.(B d — > <pK s ). We find that it is possible to satisfy 
all experimental constraints including the recent B d — > 4>K S data shown in Table IIHI in 
the framework we are considering. Furthermore, there are strong correlations between 



-4, 



CP 



, S^k and Cqk- As the accuracy of the experimental data improve, we can use 



these correlations to (in)validate the choices that we make in our model. 

Large (IrSr mixing can lead to a nonzero 9 d which depends on \5§^\ as explained in 



Section IIVBI We therefore include this new phase and perform a scan over 



32 h 



arg(5 : 



RR\ 
32 ) 



arg(/i), |5§2 | and arg(<5^ ). We show the points that satisfy all experimental constraints and 



tt The significance of the discrepancy between Sb—ts and Sb-*c is bigger, currently at about 3.5 a, where Sb^s 
and Sb^c are the averages over all measured b — > s (penguin) and b — > c modes, respectively. 
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FIG. 9: A c d ^ s7 , SfiK and C^k f° r points that satisfy all experimental constraints (within 2a), 
resulting from a scan over d^ and arg(/x), for small d^SR mixing (negligible 8d). All other 
parameters are fixed as shown in Tabled 




FIG. 10: -A^p"^" 7 , S^k and C^jf for points that satisfy all experimental constraints (within 2a), 
resulting from a scan over d^' RR and arg(/x), for large ci/js/j mixing (nonzero 9d)- All other 
parameters are fixed as shown in Table QJ 

the resulting A^p^ Xsl , and C^ K in Fig. [TUl We again see from Fig. HUlthat A^jp* Xa ^ , 
S^x and C^k are strongly correlated, although the effect of the new phase in B^B^ mixing 
allows new regions of parameter space compared to the small mixing case shown in Fig. El 
Even in the case of large mixing we find that it is possible to satisfy all experimental data 
including the S^k and C^k- One feature that we find in either large or small mixing case is 
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that sign (C^k) is positively correlated with sign (A^ d F ^ Xs ' y ). Thus, further data could shed 
light on the validity of the choices that we make in our model. 

VI. CONCLUSIONS 

A super symmetric U(2) theory has the potential to explain the gauge hierarchy and flavor 
problems in the SM. We assumed an effective SUSY mass spectrum just above the weak 
scale, the only relatively light scalars being the right handed stop and sbottom (weak scale 
masses) . We analyzed what such a hypothesis would imply for K and B meson observables by 
including all the dominant contributions that can interfere in a certain observable. Although 
for definiteness we considered a U(2) framework, our conclusions hold for any theory with a 
similar SUSY mass spectrum and structure of the squark mass matrix. 

The CP violation parameter in Kaon mixing, ex, can impose constraints on the MFV 
parameter space of our model, as we showed in Fig. |21 while the gluino contribution to ex 
is negligible. There is sufficient room to accommodate the MFV contributions to e#, given 
the present uncertainty in the lattice computation of the Bag parameter Bk- 

We find that B^Bd mixing and a^K s (sin 2/3) can impose constraints on the supersymmetric 
U(2) theory. In addition to the MFV contribution, if (IrSr mixing is large, the gluino 
contributions to BdBd mixing can be significant leading to a strong constraint on the 32 
entry of the RR squark mass matrix, 5§^, as shown in Fig. 0] Furthermore, in this case, 
there is a new phase in the BdBd mixing amplitude coming from the SUSY sector. However, 
if cIrSr mixing is small, the constraint on 5§f i from BdBd mixing is weak. 

B S B S mixing most sensitively depends on 5§^ in the SUSY U(2) theory. If 8^ is 
unconstrained by BdBd mixing (small cIrSr mixing), we showed that Am# s can be increased 
to quite large values (up to about 40 ps" 1 ), cf. Fig. |5] . The current and upcoming 
experiments can reach sensitivities required to see the SM prediction for Am^. Seeing 
a higher value, or not seeing a signal at all, might hint at some new physics of the type we 
are considering. We also presented expectations for the B s dilepton asymmetry, Af h s , which 
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can constrain 5§^. 

The experimental data on B.K.(Bd — > Xgj) imposes a constraint on the SUSY theory. 
While satisfying this constraint, we showed that an enhancement in CP Violation in — > 
X s 7 is possible. Since in the SM, A^ff* is predicted to be less than 1%, if a much 
larger value is measured, it would clearly point to new physics. In Fig. |7] we presented the 
expectations for A^, d F ^ Xal and B.K.(Bd — > X s g) while varying the magnitude and phase of 
5f 2 L . We also presented expectations for B.R.(B d — ► X s £ + £~) in Fig. |H| 

The present experimental data on the CP violation in Bd — > 01f a has about a 2 cr deviation 
from the SM prediction, and it will be very interesting to see if this would persist with more 
data. We showed that such a deviation can be accommodated in the framework we are 
considering, both for large or small (IrSr mixing. We showed, in Fig. El that A^f F ^ Xa ' 1 can 
be enhanced significantly while satisfying all other experimental bounds including the present 
data on S^k and C^k- In Figs. ElandHUJ for small and large (IrSr mixing respectively, we 
see strong correlations between A^,jT Xsl , S^k an d C$k- Comparing these with upcoming 
data with improved precision could shed light on the validity of the choices that we make in 
our model. 

We conclude by remarking that the prospects are exciting for discovering SUSY in B- 
meson processes at current and upcoming colliders. Here, we showed this for a SUSY U(2) 
model. To unambiguously establish that it is a SUSY U(2) theory, and to determine the 
various SUSY breaking parameters, will require looking at a broad range of observables. 
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APPENDIX A: MIXING ANGLES 



The charged SU(2) Majorana gauginos Wi, W 2 can be combined to form the Dirac spinor 

W^U W \ ) , (Al) 

where W^ 1 = W± a ± iW 2a - The up and down type Higgsinos can be combined to form the 
Dirac spinor 

The chargino mass terms can then be written as 



£D-(W+ H + )(M x P L + M] ( P R )l f , + j, (A3) 
where 

/ M 2 V2sm3m w \ 

M * = \m r )■ (A4) 

\ V 2 cos p mw A* / 

We can go to the chargino mass eigen basis (xi X2) by making the rotations 

MI+) =(Ci ' fl) M*D , (A5> 

with the rotation matrices C^ tR given as 

/ cos6 a - sin 9 a e~ ila \ ( e iria \ 
C a = . , (A6) 

\smB a e 1 ^ cos6 n J V e lpa 

where the mixing angles and phases are 



tan 20 



7 L = - arg (M 2 + fi cot (3) (A7) 
mw sin (3\M 2 + fJ, cot (3\ 



L 



tan2# R 



M| + |/i| 2 + 2 cos 2/? 
j R = - arg (M 2 + /i* tan /3) 

8mvi/ cos/9|M 2 + /z*tan/?| 



M|- |/i| 2 -2m^cos2/? 
7^ = arg (c R (M 2 c L + v^m^ shi/SsLe* 71 -) + Si J e" i7 ' R (v / 2m H / cos j3c L + /is^e* 

p fl = arg (c R (-y/2m w cos (3s L e~ nL - fiC L ) + SKe nil (-M 2 SLe" m + y/2m w sin/5ci) 
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with < 9 a < n / 2 so that M% > M% 2 , and, sl,r = sin^jj and cl,r = cos 9 l,r- 
The sbottom mass terms are given as, cf. Eq. (JSJ 

£ = f + + ^ w»»n ,5^ (A8) 

Vud-Aj - /i*tan/3 m 6 m\ R + m 2 , + J \b R J 

where the A^ are the D-term contributions given as 

A d = (T 3 - Qem sin 2 9 W ) cos2/?m| . (A9) 

This is diagonalized by the rotation 

fb L \ ( cos9~ b -sin^e-^\ /6A 
UjJ Uin^e^ cosflg / \b 2 ) 

where the mixing angle and phase are given by 

2|t>dv4f ) — /i* tan (3 rrib\ 



tan 26* 



6 ~ (rnj LL + Ai; 
7j, = arg(f d A b - /i* tan/5 m 6 ) . (All) 



(rnz LL + A d L )-{ml R + A d R ) 



We have similar equations for stop mixing with obvious changes, in addition to the off 
diagonal term now being given as: (v u A t — //* cot (3 m t ), and the stop mixing matrix denoted 
as Cf. In our framework, owing to the smallness of the off diagonal RL mixing term compared 
to rn\ LL ~ m 2 ,, we have small stop and sbottom mixing. Furthermore, the sbottom mixing 
angle is negligibly small and we neglect its mixing effects. We thus have b\ ~ b^ and bi ~ 6r. 
The stop mixing angle, however, is not as small and so we include its effects. 

To compute the interaction vertices in the SuperKM basis, one could diagonalize the 
6x6 squark mass matrix. Since the off-diagonal entries in our case are small, we perform 
an approximate leading order diagonalization of the mass matrices shown in Eq. (jSJ . 

Focusing first on the &r§l mixing, 

c-isi h) (^^;^t ^* )('*). ( A12) 



b R 
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This is diagonalized by the rotation 

/s L \ ( cos9§t - sin B^e^ \ / q 2 
\b R ) ~ [smO^e^ cos^ 2 L ) \q 3 , 



-(<U.) * V (A13) 



The mixing angle and phase are given by 
tan2#£ L 



2\v d Ae\ 



:;2 ' [m\ LL + ml + A£) - (m? fl + m\ + A&) ' 



7 3 T = arg(^e) . (A14) 



The b^SR mixing is given similarly. We have the mass terms 

V emj mf fl + + A£ / \ 6 fl / 

which is diagonalized by the rotation 

/s R \ I cos 6™ -sin^e-^X /<? 2 \ 
UjJ \sin0$VT£ H cos^ 2 R /Us/ 

- C) • 

The mixing angle and phase are given by 

tai1 32 ' + ml + Aj) - (m? fl + + A|) ' 

7^ = arg(em2) . (A17) 



The <1lSl and d,RSR mixing terms are 

/ m\ + m\ + A d L R ie'ml \ ( d L ,n \ 

£D-(d* LR r LR )l ' • ( A18 ) 

V -ie'ml ml + m 2 s + e 2 ml + A d LR ) LLRR \s L , R ) 

The matrices that diagonalizes these, Cj L g L and C^ rSr are given analogous to Eq. flHBJ), the 
angle and phase (#f 2 L , I12) an d ($12^ 7r2 R ) given analogous to Eq. (jA17|) . and we will not 
write them down explicitly. The diagonal entries are split only by 0(e 2 ), and therefore this 
mixing is maximal in general. However, if e'mj <C e 2 m| this mixing can be small. 
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APPENDIX B: LOOP FUNCTIONS 



The Bd — > X s 7 loop functions are given by 

, LL (x(7-5x-8x 2 ) x 2 (3x-2) 



F 7 LL (x) 
F 8 LL (x) 
F 7 RL (x) 
F 8 RL (x) 
F 7 LL (x) 
F* LL (x) 



+ 



3Q(x-l) 3 6(a;-l) 4 

'x(2 + 5x- x 2 ) 3x 2 

12(x - l) 3 Q(x- l) 4 

5 - 7a; a;(3a; - 2) , \ 

mxl , 



In a; 



In a; 



+ 



6{x-l) 2 3(a;-l) 3 
1 + x x 



2(x-\) 2 {x-iy 

x(3 — 5x) x(3x — 2) 



+ 



12(x-l) 2 6(a;-l) 3 
a;(3 — x) x 



In a; J , 
In a; I , 



F,(x) = 
F~ g (x) = 



A(x-l) 2 2(x-lf 
x 2 — 1 — 2a: In x 



In a; 



2(x - l) 3 
5a; 2 -18a; + 13 + (9 - x) In x 
3{x- l) 3 



The SM box functions are given by 

4a; — 1 la; 2 + x 3 3x 3 In x 

So{x) = 



S (x u x 2 ) = | 



4(1 - a;) 2 2(1 - a;) 3 ' 
X\X 2 f x\ — 8a;! +4 x\ — 8a; 2 + 4 



{xi -x 2 )(x 1 - l) 2 

3 — ] 

-1)1 



lnxi + 



(x 2 - Xi)(x 2 - 1)' 



\nx 2 



( Xl -l)(x 2 -l) 
The charged-Higgs and chargino box functions are given by 



In r a + 



(Bl) 



(B2) 



(B3) 



■In 773 



{rp - r a ) (si - r a ) (sj - r a ) (r a - rp) (s, - rp) (sj - r p ) 

s 2 s 2 
+ - — 1 — — ^ In Si + — — ^ In sj , 



(r a - Si)(rp - Si)(sj - s^ 



Y 2 {j' w rp, Si, Sj) — y/ SiSj 



In r a + 



(r a - Sj )(rp - Sj)( Si - Sj ) 

r f3 



(rp - r a ) (si - r a ) (sj - r a ) (r a - rp) (si - rp) (sj - rp) 



In rp 



In Si + 



r a - Si)(rp - Si)(sj - Si) (r a - Sj)(rp - Sj)(si - Sj) 



In s. 



(B4) 
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from which various limiting cases can be obtained. 



The gluino box integrals are given as 

d 4 p 



h = h(M? g ,Ml,ml,ml R ) = J 



h = h{MlMlml R ,ml Ri 



(2vr) 4 (p 2 - Mf)(p 2 - Mf)(p 2 - m 2 b J(p 2 - m? ) 
d 4 p 



P 



(2vr) 4 (p 2 - M?)(p 2 - M?)(p 2 - mf )(p 2 - m\ ) 
which we evaluate numerically using LoopTools in Mathematica. 
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